Anomalous superfluidity in 2 + 1 dimensional two-color lattice QCD 
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We study thermodynamics of strongly coupled lattice QCD with two colors of staggered fermions 
in (2 + 1) dimensions. The partition function of this model can be written elegantly as a statistical 
mechanics of dimers and baryonloops. The model is invariant under an 50(3) x (7(1) symmetry. 
At low temperatures we find evidence for superfluidity in the U(l) symmetry sector while the 
5*0(3) symmetry remains unbroken. The finite temperature phase transition appears to belong to 
the Kosterlitz-Thouless universality class, but the superfluid density jump p s (T c ) at the critical 
temperature T c is anomalously higher than the normal value of 2T c /ir. We show that by adding 
a small SO(3) symmetry breaking term to the model, the superfluid density jump returns to its 
normal value implying that the extra symmetry causes anomalous superfluid behavior. Our results 
may be of interest to researchers studying superfluidity in spin-1 systems. 

PACS numbers: 



Two dimensional field theories are dominated by 
infrared fluctuations which lead to interesting non- 
perturbative phenomena. For example consider the 
O(N) non-linear sigma model in two dimensions : When 
N > 2 the theory is known to be asymptotically free 
and develops a non-perturbative mass gap, while the 
N = 2 theory has both a critical (superfluid) phase and 
a massive (normal) phase, separated by the well known 
Kosterlitz-Thouless phase transition The transition 
is topological in nature and is driven due to the unbinding 
of topological defects (vortices). One of the striking the- 
oretical predictions is that the superfluid density p s (T c ), 
at the critical temperature T c , is known to jump from 
2T c /ir to zero and has been observed in superfluid 
helium Q and recently in spinless atomic BEC Q. 

Recently theories with an 5*0(3) x U(l) symmetry have 
become interesting due to the possibility of creating Bose 
condensates of atoms with non-zero spin, such as 23 N a, 
39 K and 87 Rb, in quasi- two dimensional atomic traps li- 
lt is interesting to understand if these exotic spin systems 
can show new features. In a recent work it was argued 
that due to the presence of half vortices, superfluidity 
arising due to the condensation of spin-1 particles will be 
anomalous and in particular the superfluid density jump 
at T c will be 8T c /tt instead of 2T c /ir |f|. Some evidence 
for this conjecture was presented using Monte Carlo cal- 
culations on lattices as large as 40 2 . The Hamiltonian 
used to discuss the spin-1 spinor condensation was first 
introduced in 0, where it was argued that the ground 
state can either be in a "polar" state or a "ferromag- 
netic" state. The fundamental difference between these 
two ground state manifolds leads to a different vortex 
structure such that superfluidity without spin order ap- 
pears possible only in the polar phase 0, S B ■ 

In this article we show that anomalous superfluidity 
also emerges naturally in a 2 + 1 dimensional strongly 
coupled lattice gauge theory where fermions interact with 
an SU (2) gauge field. As we will argue below, this model 
also has an SO(3) x U(l) symmetry and hence may be 



useful as an alternative model to study superfluidity in 
spin-1 systems. Fortunately, the partition function can 
be rewritten elegantly as a statistical mechanics of dimers 
and baryonloops making our model accessible to a vari- 
ety of researchers. In this representation our model can 
be studied efficiently with ease on lar ge la ttices, thanks 
to recent developments in algorithms [10(. Here we re- 
port on results obtained from studies on lattices as large 
as 1024 2 . We find clear evidence for anomalous super- 
fluidity: the jump in p s (T c ) is not close to its normal 
value of 2T c /ir. On the other hand, it is neither close to 
8T c /ir as conjectured in Q making our result even more 
interesting. Two-color lattice QCD (2CLQCD) has been 
studied in 3 + 1 dimensions extensively as a toy model for 
the physics of strong interactions at finite temperatures 
and densities and continues to be interesting even today 
[ll| . In the strong coupling limit it was studied recently 
using mean field theory [12j and Monte Carlo methods 
|l3| . The theory has remained relatively unexplored in 
(2 + 1) dimensions (however see 0]). 

The action of 2CLQCD that we study here is given by 
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where x represents a site on a cubic lattice with coor- 
dinates (x t ,xi,x 2 ), = (x(x),~X tr ( x ) T 2), and ^ x r = 
(x tr ( x )i x( x ) T 2) are four component Grassmann fields as- 
sociated with even and odd lattice sites respectively. In 
our notation a are Pauli matrices that mix x and x tr 
present in '5 and \& while r are Pauli matrices that act 
on the color space. The Grassmann valued quark fields 
xi x ) and x( x ) represent row and column vectors with 2 
color components. The gauge fields U a . x are elements of 
SU{2) group and live on the links between x and x + a 
where a = t, 1, 2. The factors r)t, x — 1, r]i_ x = (— l) Xt and 
i]2, x = {— l) xt+xi are the staggered fermion phase factors. 
We fix a\ = oi2 = 1 but vary a t . By choosing a slab ge- 
ometry of the type L t x L 2 (periodic in all directions) 
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FIG. 1: An illustration of a dimer-baryonloop configuration 
in (1 + 1) dimensions. 



we can study thermodynamics in the L — > oo limit at a 
fixed L t . Our results were obtained with L t = 4. The 
parameter S = l/(a 4 ) 2 then controls the temperature. 

The action JTJ is invariant under <3> — * V$, ¥ — > fV^ 
where V = exp(i</> • a + iff) € t/(2). Since V — 1 when 
= (0, 0, 7r) and 9 = 7r the real symmetry of the model 
is 50(3) x U(l), similar to the symmetry of the model 
studied in [6j- In two spatial dimensions, at zero temper- 
atures, the above model is expected to develop a "chi- 

f T 

ral condensate" which means {^ tr a\^) = (^oi^ ) = 
(XX) 0. Clearly this breaks the 5*0(3) x U(l) symme- 
try to 0(2). Note that the order parameter is invariant 
under V = exp(i03<73)(l, <7i) which is simply isomorphic 
to 0(2). Thus, the Goldstone boson manifold is the same 
as that of the polar phase discussed in [6J, |9j . This con- 
nection should make our work interesting to researchers 
studying spin-1 spinor condensates. 

It is possible to integrate out the gauge fields exactly 
and write the partition function as a statistical mechan- 
ics of dimers with constraints [l3l Hj| . A dimer is a bond 
variable connecting two neighboring sites x and x + a. 
There are two types of dimers possible in our model: 
an oriented dimer (with an arrow) and a dimer with- 
out orientation. Oriented dimers will be represented by 
b a {x) = 1,-1, while non-oriented dimers will be repre- 
sented by d a {x) = 1,2. If a particular type of dimer is ab- 
sent then the corresponding b a (x) or d a (x) is set to zero. 
Clearly d a (x) — d- a {x + a) and b a (x) — —b- a (x + a). 
The microscopic model also introduces two constraints: 
(1) the sum of d a (x) + \b a (x)\ over all the bonds a con- 
nected to a site x must be 2, (2) every site must have 
one incoming and one outgoing oriented dimer or no 
oriented dimers. Thus, each configuration consists of 



three types of connected objects: (a) self avoiding ori- 
ented loops made with oriented dimers (called baryon- 
loops), (b) self avoiding non-oriented loops made with 
non-oriented dimers, and (c) double dimer bonds. For 
illustration we show an example of a (1 + 1) dimensional 
dimer-baryonloop configuration in figure ^ However, 
note that our work is done in (2 + 1) dimensions. 

The Boltzmann weight of each configuration is given 
by 8 n where n is the total number of temporal dimers 
(oriented or otherwise) . We focus on two types of winding 
number susceptibilities: 



(2a) 
(2b) 



constructed from the conserved currents Ja^~\x) = 
cr x [\b a (x)\ + d a {x) — 2], and Ja°^ = b a (x) which arise 
due to the U(l) and 5*0(3) symmetries discussed above. 
Here, a x — 1 on even sites and —1 on odd sites. 

In the superfluid phase we expect Xw ^ 7^ for large L. 
In figure |21 we plot Xw as a function of L for different 
values of 6. For 5 < 0.6, Xw^ ^ s independent of L in 
the range 32 < L < 1024, suggesting that we are in the 
superfluid phase. When S > 0.8 we see that Xw begins 
to decrease for large L suggesting that superfluidity is 
lost. For 5 = 0.7 our data fits well to the form 



X 



u(l) 
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1 



0.5 



log(L) + B 



(3) 



expected to hold at the KT transition. We find A = 
1.59(1) and B = 2.7(7), with a x 2 /DOF = 0.5. Since 
= Ps{T)/T one would have expected A to 
be close to either 2/-7T for normal superfluidity or 8/tt 
for superfluidity due to half vortices |6j. Our data is 
inconsistent with both these expectations. 

In two dimensions there is no spontaneous symmetry 
breaking which means we should expect Xw to go to zero 
for large L. Of course U(l) symmetry is special which 
allows for a non-zero x^ 1 ^, but Xw^ must vanish for 
large L. If £ is the correlation length in the 50(3) chan- 
nel, we expect Xw^ is exponentially small for L » 
However, for small values of S our results are obtained in 
the region where L <C £ since £ is expected to be asymp- 
totically large. Here the L dependence of Xw^ can be 
found by solving a one-loop renormalization group equa- 
17] which gives 



tion 



X 



~W =C-Dlog(L). 



(4) 



Our results for 6 < 0.7 are consistent with this behavior 
while at 8 = 0.7 it just begins to break down suggesting 
that £ ~ 1024 there. The fits are as follows: 
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FIG. 2: Plot of X™ versus L for different values of 8. At 
8 = 0.1,0.5 and 0.6 we find x l W = 2.335(3), 1.963(3) and 
1.838(3) fits the data well (solid lines). At 8 = 0.7 eq.J3J with 
A = 1.59(1) and B = 2.7(7) fits better. (dashed line). 
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FIG. 3: Plot of Xin^ versus L for different values of 8. The 
solid lines are fits to eq.l^J discussed in the text. 
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For larger values of 8, Xw^ begins to decrease faster. 
These results are plotted in figure EH 

Based on the above analysis we conclude that in our 
model the U(X) symmetry leads to anomalous superfluid- 
ity at low values of S while the SO(3) symmetry remains 
unbroken. We estimate that 6 C is roughly between 0.6 
and 0.8. An obvious criticism to our conclusions is the 
following : Perhaps we have not reached the thermody- 
namic limit for S ~ 0.6 since in this region L <C £■ Is it 




FIG. 



4: Plot of xl W versus 8 at L = 256 for W = 1 (tri- 
angles) and W — 0.9 (squares). The solid lines just connect 
the data points to guide the eye. The inset shows the finite 
size scaling plot similar to figure H at W = 0.9. The values 
2u/n, v = 1, 2, 3, 4 are shown as dashed lines for reference. 



L ~ £ and then begins to drop to a different value |18|? 
If true, this would imply that 5 C could be much smaller 
than our estimate so that the jump in Xw^ at S c could 
indeed be 8/71". 

If the above criticism is true, one would expect that 
Xw is very sensitive to £ in the region S ~ 0.6. We 
should be able to study this sensitivity by changing £. 
We can do this by changing the Boltzmann weight of each 
dimer-baryonloop configuration to 5 n W m where m is the 
number of oriented dimers in the configuration. When 
W = 1 we obtain the original model, but when W < 1, 
£ is considerably shortened since the SO(3) symmetry of 
the original model is explicitly broken to a U(l) subgroup 
[l9|. Changing W does not affect the superfluid U(l) 



symmetry. The two currents J«° and continue to 
be conserved due to the remnant U(l) x U(l) symmetry 
of the modified model. 

The W — model was studied in where it was 
shown that Xw"^ shows the normal 2/7T jump at the KT 
transition. When W = there are no oriented dimers 
which means that x^ 3 '' = exactly. In this work we 
have studied the modified model for W = 0.9 and we 
again find that Xw^ = for the lattice sizes we have 
studied showing again that £ has shortened considerably. 
In figure 0] we plot Xw as a function of 8 for L = 256. 
For comparison we plot the values for both W = 0.9 and 
1. The inset of figure0]shows the dependence oixw^ 011 
L at 8 = 1.5 and 1.6 for W = 0.9. The data at 1.5 fits to 
the form given in eq.JSJ, with A = 0.771(4), B = 2.2(4) 
with x 2 /DOF = 0.7 while at 1.6 the fit is poor. This 
suggests 1.5 < 5 C < 1.6 and the jump in Xw at 8 C is 



not possible that Xw remains almost a constant until now very close to the normal value of 2 /it. 
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FIG. 5: Plot of <^ 1 (extrapolated to the infinite volume) 
versus S for W = 1. 

Clearly, the breaking of the SO (3) symmetry has a 
significant effect on the critical point and the jump in 
the superfluid density at the critical point returns to its 
normal value. However, changing W from 1 to 0.9 has 
negligible effect on Xw in the region 8 ~ 0.6. This 
weakens significantly the criticism raised above that our 
lattices are far from the thermodynamic limit for 8 ~ 0.6. 

We have also measured the correlation length £ us- 
ing the correlator (x x Xx X v Xv) where x and y are two 
lattice sites. We look for an exponential decay along 
a spatial direction while projecting to zero momentum 
(or frequency) in the other directions. Figure [5] shows 
the behavior of C _1 extracted from our data after re- 
moving all finite size effects. We see that C appears to 
be diverging rapidly. However, note that the correla- 
tor we measure couples to fluctuations which have both 
SO(3) as well as U(l) quantum numbers. Hence it is 
possible that ( grows asymptotically large as 8 is low- 
ered but does not diverge at the critical point. Clearly 
our data cannot be used to confirm this. However, in 
order learn something, here we fit the data to the form 
= a exp(— (3/^/8 — 8 C ), expected at a KT transition. 
We find a = 47(15), (3 = 4.59(3) and S c ~ 0.595(20) with 
a x 2 /DOF — 0.85. Interestingly this estimate of 8 C is 
not very different from our previous estimates. 

Before we end we would like to note that it is possible 
to modify our model so that the jump at T c is 8T c /n as 
predicted in 0. For example if the Boltzmann weight 
of each dimer-baryonloop configuration is chosen to be 
8 n (WD) m , where m is the number of d a (x) = 2 bonds 
in the configuration, then it is easy to argue that the 
50(3) x U(l) symmetry remains unaffected. For large 
Wd all loops will be suppressed and the model should 
behave like the close-packed dimer model studied in • 
The only difference is that the current will now 

be twice the value in the earlier study. Since the earlier 



study showed a normal jump in the superfluid density, 
our model with large Wd should show four times the nor- 
mal jump. We have confirmed this for Wd = 3.5. This 
is not surprising since increasing Wd makes the SO(3) 
symmetry more disordered and thus makes half-vortex 
excitations easy. It would be interesting to study how 
the superfluid density jump evolves with Wd and under- 
stand if an interplay between half-vortices and normal 
vortices can lead to the anomalous superfluid behavior 
observed here. 
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